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ABSTRACT 

Using all available multi-component radio pulse profiles for pulsars with medium to 
long periods and good polarisation data, we have constructed a two-dimensional image 
of the mean radio beamshape. This shows a peak near the centre of the beam but is 
otherwise relatively uniform with only mild enhancements in a few regions. This result 
supports the patchy-beam model for emission beams in which the mean beam shape 
represents the properties of the emission mechanism and observed pulse components 
result from emission sources distributed randomly across the beam. 
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1 INTRODUCTION 

Pulsars are generally believed to be rotating neutron stars 
in which the observed pulses result from one or more emis- 
sion beams sweeping across the E arth as the star rotate s. 



Obse rvations of radio polarisation ( Radhakrishnan & Cooke 



196E ) led to the magnetic-pole model in which the emission 



beam was centred on the magnetic axis of a predominantly 
dipole magnetic field. On the assumption that the emission 
is directed radially, the observed pulse profile reflects the 
variations in emission intensity along a line of constant ro- 
tational katrtud£ ; _Alth£>u^ 

known ( Taylor. Manchester Lyne 1993} Lyne et, aj 200f 1 



in the plane of curvature of the magnetic field and hence the 
model naturally explained the smooth sweep of polarisation 
position angle seen in many pulsars, as well as the common 
occurrence of double-peaked pulse profiles. 

Backer (1976) first discussed the more-or-less central 
pulse component seen in many pulsars and suggested that 
it resulted from an axial beam. Rankin (1983) made the 
distinction between this central or 'core' emission and the 
outer or 'conal' emission, and showed that these two compo- 
nents of the profile had rather different properties. At least 
in short-period pulsars, cor e emission ha s a steeper spec- 



the sha pe of pulsar radio beams remains controversial. 



trum than conal emission ( |Rankin 1983] ; [Lyne fc Manch- 



ester 1988) and the fluctuati on properties of the two types of 



There are two main areas of uncertainty. One concerns 
the outline shape of the radio beam. Originally assumed to 
be circular, some investigations (e.g. Narayan & Vivekanand 
1983) argued for an elliptical beam extend in the latitude 
direction with a large axial ratio. In contrast, Biggs (1990) 
suggest ed that the beam was slightly compressed in the lat — 
itudinaLdixeclkirj Moat other inaestigatknis ( p g U mp Rr 



emission are very different (Rankin 198£). This model was 



subsequently extended to have two or more coaxial conal 
emission zones, either to account for multiple-com ponent 
profiles (jGlil fc Krawczyk 199?]; |Qiao fc Lin 1998j) or for 



the appearance of c omponents at different apparent r adii 
from the conal axis (Rankin 1993a Kramer et al. 1994 



Mi- 



tra & Deshpande 1999). These latter analyses depend on a 



Manchester 1988; Bjornsson 1998; Gil & Han 1996) have 
concluded that the emission beam is essentially circular, and 
we will assume this in the present investigation. 

The second area of uncertainty concerns the form of the 
beam pattern. Early observations showed that there are two 
or more pulse components in many pulsars (e.g. Lyne, Smith 
& Graham 1971; Manchester 1971). Ruderman & Suther- 
land (1975) presented a detailed model for pulsar radio emis- 
sion in which the beam had the shape of a hollow cone, more 
intense around the periphery, corresponding to the last open 
field lines emanating from the polar cap region. The radio 
emission was attributed to curvature radiation by positrons 
moving these field lines. Such radiation is linearly polarised 



knowledge of the angle between the beam and rotation axes 
(a), usually computed from the observed width of a 'core' 
component. Unfortunately, the determination of this angle is 
very model dependent and this results in large uncertainties 
in the derived radii. 

Based on a large sample of pulse shape and polarisa- 
tion data, Lyne & Manchester (1988) found that pulse com- 
ponents in multiple-component profiles were not symmetri- 
cally located within the beam boundary. Furthermore, they 
were generally of very different intensities and, in some cases, 
missing altogether. Their results were consistent with com- 
ponents being randomly located within the beam bound- 
ary, leading to the 'patchy-beam' model. Gould (1994) and 
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Table 1. Normalized impact angles and the frequency of profile samples 
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Gould, Lyne & Smith (2000) confirmed these conclusions 
with an even larger data set. 

Manchester (1995) suggested an interpretation of these 
results in which the observed pulse profile is the product of 
a 'window function', which is a function of pulse period and 
radio frequency but common to all pulsars, and a 'source 
function' which determines the strength of the emission at a 
given point within the beam and is different for every pulsar. 
For example, the source function may be determined by the 
density or energy of the plasma beams along different field 
lines, whereas the window function is determined by the 
properties of the emission mechanism. 

In this paper, we determine the two-dimensional shape 
of the window function at frequencies around 1 GHz by av- 
eraging the observed pulse profiles of multi-component pul- 
sars having reliable polarisation data. This analysis depends 
only on the observed width of the pulse profile and the 'im- 
pact parameter', that is the minimum angle between the 
line of sight and the beam axis, usually given the symbol 
P. In contrast to the beam inclination angle a, the impact 
parameter is normally well determined from the m aximum 
rate of change of p osition angle at the profile centre ( Lyne & 
Manchester 198S ) . We assume that emission is seen across 
the entire polar cap and choose only pulsars with profiles 
where this appears to be the case. In Section 2 we de- 
scribe the data set and our methods of determining the 
two-dimensional beam pattern. Results are presented and 
discussed in Section 3 and Section 4 summarises our conclu- 



2 THE DATA SET AND ANALYSIS METHOD 

In order to reliably determine the two-dimensional beam 
shape we choose pulsars which have pulse profiles with two 



or more pulse components, good signal/noise ratio and good 
quality polarisation data. 

For most pulsars with two or more pulse components, 
we see emission from right across the polar cap. The polari- 
sation information allows a measurement of the normalised 
impact angle, /3 n , that is, the impact parameter (3 expressed 
as a fraction of the beam radius. This parameter is com- 
puted assuming a — 90°, but it is not very sensitive to the 
actual value of a (Lyne & Manchester 1988). The polari- 
sation data also allow a check on the assumption that we 
are seeing emission right to both edges of the beam. Only 
those pulsars for which the maximum rate of change of po- 
sition angle is close to the pulse centre were included in the 
sample. Single component pulsars are not included in the 
sample as it is often difficult to reliably determine if they 
are predominantly core emission (low /3 n ) or from the outer 
edge of the beam (high /3 n ). 

Millisecond pulsars were also excluded from the sample 
as they generally have very wide and complex profiles and 
polarisation variations which do not fit the simple rotating- 
vector model (e.g. Navarro et al. 1997, Stairs et al. 1999). 
It is therefore difficult or impossible to obtain reliable val- 
ues of j3n for these pulsars. Very short-period pulsars also 
commonly have wide or interpulse profiles and are excluded 
from the sample for similar reasons. 

Pulse profiles in a given pulsar usually evolve with fre- 
quency, with some components becoming stronger and other 
becoming weaker or even disappearing at higher or lower fre- 
quencies. At low frequencies, the core component tends to 
dominate the profile, whereas at high frequencies the conal 
emission is generally dominant. To ensure relative unifor- 
mity in the data set, we restricted our attention to pro- 
files at frequencies between 600 and 1600 MHz. The results 
we obtain therefore represent pulsar emission at frequencies 
around 1 GHz. 

We checked the profiles available on the pulsar profile 
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Normalised Pulse width 



Figure 1. Beam components for four pulsars derived using the 
procedures described in Section 2. The pulse profiles arc given in 
the lower panels where the dotted lines represent the interpolated 
profile used in the projection on to the beam pattern. 



database ^ of European Pulsar Network maintained at Max- 
Planck-Institut fur Radioastronomie (Lorimer et al. 1998) 



which contains pulsar profiles from more than 50 papers, 
including the large datasets from Gould & Lyne (1998), 
Manchester, Han & Qiao (1998) and Hoensbroech & Xilouris 
(1997). Normalised impact parameters were obtained from 
Lyne & Manchester (1988), Rankin (1993), Gould (1994) 
and Manchester et al. (1998). The final data set of 87 pul- 
sars which satisfy the above criteria is listed in Table 1. 

To determine the average beamshape, we first defined 
the duration of each pulse profile by taking the points at 




-0.5 0.5 

Normalised beam radius in longitude 



Figure 2. The beam pattern obtained by adding data for all 
pulsars in the sample represented as a greyscale. This pattern 
has not been normalised for the nonuniform distribution of j3 n . 



3 



0.2 0.4 0.6 0.8 1 

Normalised impact angle /? n 

Figure 3. Distribution of normalised impact paramters, j3 n , for 
pulsars in the sample. The dotted line is the sum of the gaussian- 
broadened values used in the projection (see text). 



which the signal rose above 3a, where a is the rms devia- 
tion of the off-pulse noise. Next, the mean intensity of the 
observed emission within the pulse was normalised to unity. 
The normalised impact parameter, (3 n , gives the offset of 
the locus of line of sight across the polar cap in the latitude 
direction as a fraction of the beam radius. We assumed a 
circular beam of unity radius and a straight trajectory of 
the emission point across the polar cap. The beamshape is 
assumed to be symmetric about the equator, i.e., the sign 
of /3 n is ignored. 

The profile intensity was interpolated on to an x — y ar- 
ray representing a semi-circular beam of unity radius along 
a line at y = f3 n between x values of ± \/ 1 — [3^ using a five - 
point polynomial interpolation routine (Press et al. 1992). 



To allow for uncertainties in the value of /3 n , possible non- 
linearities in the beam trajectory across the polar cap, and 
to represent the finite width of a subpulse beam, the pulse 
profile is broadened in latitude (y) assuming a guassian form, 
exp[— (/? — /3 n ) 2 /0.1], between A(3 = ±0.3. Examples of com- 
puted beam components for two pulsars are shown in Fig. 
1. 



See http:/ /www. mpifr-bonn.mpg.de/div/pulsar/data/ 
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Figure 4. Average shape of the pulsar radio beam for /3 n < 0.8. The contours are at multiples of 0.1 of the peak value near the beam 
center. 



3 RESULTS AND DISCUSSION 

Integrating all pulse profiles in the sample gives the beam 
shape shown in Fig. 2. This figure represents the probabilty 
of observing beam components at each location. However, 
it is strongly affected by the sample distribution of f3 n . As 
shown in Fig. 3, this has a strong peak at /3 n ~ 0.25, indi- 
cating that multi-component pulse profiles are most likely to 
be detected at about this impact parameter. The deficit at 
f3 n ~ is an observational selection effect which results from 
the smearing of rapid position angle changes near the profile 
centre due to finite subpulse beamwidths and instrumental 
broadening of the profile. Few multi-component profiles are 
observed at high /3„. 

The beam distribution shown in Fig. 2 was normalised 
to correct for the non-uniform distribution of /3 n using the 
gaussian-smoothed distribution shown in Fig. 3, giving the 
final average beam shape shown in Fig. 4. Since there are 
only one or two pulsars in bins with j3 n > 0.8 (Fig. 3), there 
is large statistical uncertainty in the average profile shape 
for these bins. Therefore we do not plot this region of the 
normalised beam pattern. 

The overall impression given by Fig. 4 is that of a 
'patchy' beam. There are some systematic features, but in 
general the intensity is relatively uniform across the whole 
beam. The central or 'core' feature is significantly displaced 
toward later longitudes. The next most prominent property 
of the average beam after the core component is the pres- 
ence of a few patches of slightly enhanced emission. Most of 
these are located within an irregular and rather broad an- 
nular maximum at a normalised radius of about 0.7 which is 
somewhat stronger on the leading side of the beam. Under- 
lying this is more-or-less uniform emission over the whole of 
the emission beam. The decline at the beam edge is not very 
sharp. 



Although there is no evidence for double or multiple 
cones in Fig. 4, our present data set is dominated by two- 
and three-component pulsars. A larger sample of pulsars 
with more than three components is needed to reliably dis- 
tinguish multiple-cone models from a patchy beam models. 
However, our results suggest that, if multiple cones exist, 
they are at different radii relative to the beam radius in dif- 
ferent pulsars. They also show that the conal emission is not 
confined to a single annular region at the beam boundary. 

These results are consistent with the idea that compo- 
nents in pulse profiles are largely determined by the distri- 
bution of sources across the polar cap - the 'source function'. 
The rather smooth distribution distribution of intensity in 
the mean beam shape suggests that these source regions 
are randomly distributed for different pulsars as suggested 
by Lyne & Manchester (1988) and Manchester (1995) - 
the 'patchy-beam' model. The number of identifiable source 
components across the profile is mostly limited by the finite 
width of subpulse beams. Except in a few c ases with ver y 
high signal/noise ratio, e.g. PSR B0740-28 (|Kramer 199^ ), 



it is generally not possible to identify more than four or five 
peaks or components, and often only one or two. Some very 
wide components identi fied by gaussi an fitting, e.g. PSRs 
B2319+60, B2021+51 ( |Kramer 1994j ), are most probably 



regions of distributed emission producing overlapping sub- 
pulse beams. 

We believe Fig. 4 is a good representation of the mean 
radio beam (for f3 n < 0.8) emitted at frequencies around 
1 GHz by 'typical' pulsars, that is, the pulsars of medium 
or long period which dominate the sample of known radio 
pulsars. It represents the 'window function' in the model of 
Manchester (1995) which is determined by the effective gain 
or efficiency of the radio emission process. 
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4 CONCLUSIONS 

We have computed the average radio beamshape at frequen- 
cies about 1 GHz of pulsars with medium to long periods. 
This beamshape has a peak near its centre and a mild, broad 
and rather irregular enhancement at a normalised beam ra- 
dius of about 0.7, but is otherwise rather uniform. The de- 
cline at the beam edge is gradual. 

These results suggest that the presence and location of 
profile components are determined by a 'source function' 
which varies randomly from pulsar to pulsar. The summing 
of these randomly distributed components results in a rela- 
tively uniform average beam profile which we interpret as the 
'window function' representing the properties of the emis- 
sion process common to all longer-period pulsars. 
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